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Abstract 

We introduce models of generic rigidity percolation in two dimensions on hierarchical networks, 
and solve them exactly by means of a renormalization transformation. We then study how the 
possibility for the network to self organize in order to avoid stressed bonds may change the phase 
diagram. In contrast to what happens on random graphs and in some recent numerical studies at 
zero temperature, we do not find a true intermediate phase separating the usual rigid and floppy 
ones. 
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I. INTRODUCTION 



Consider a structure made of sites, connected by links. Each link imposes a constraint, 
by prescribing the distance between the two sites it connects: if the actual distance between 
the two sites is different from the prescription, there is an associated energy cost. Rigidity 
theory deals with properties associated with the structure's topology, which do not depend 
on the physical nature of the constraints, nor on the precise form of the energy cost. Typical 
questions asked by rigidity theory are: how many degrees of freedom are left in the system? 
Is there a macroscopic cluster of sites rigidly connected one with the others? By contrast, 
questions related to the elastic properties of the structure do depend on the specification of 
the constraints. When the number of links in the structure is increased, the phenomenology 
is as follows. For a small enough number of links, that is a small mean connectivity of the 
structure, there are many more degrees of freedom than constraints; there is no macroscopic 
rigid cluster, and many degrees of freedom are left in the system: the system is said to 
be floppy. At large mean connectivity, there are many more constraints than degrees of 
freedom; there is a macroscopic rigid cluster, and many constraints cannot be satisfied: the 
system is said to be stressed rigid. In between these two phases takes place the rigidity 
percolation transition. 

Despite its clearly mechanical origin, the problem of rigidity percolation has also attracted 
attention in the last 30 years because of its applications in understanding the properties of 
network forming glasses, like GeSe or GeAsSe alloys 

A3 

. In this case, the atomic bonds 
may be considered as constraints. However, contrary to standard rigidity percolation as 
presented in the first paragraph, angles between adjacent bonds also need to be considered 
as constraints, in addition to bond lengths. In the following, we will concentrate on central 
force rigidity (that is when only bond lengths are considered as constraints), in two dimen- 
sions. Beyond the applications to glasses, models of cross linking stiff fibers forming random 
networks were shown to fall in the same universality class as central force 2D rigidity 



this type of system has been used to model network forming proteins 

A decisive theoretical progress was made in the 90's for generic {37I 2D rigidity, with the 
introduction of combinatorial algorithms 6|, |7|], based on Laman's theorem [8|. These new 
algorithms allow for the study of much larger samples than before, as well as more precise 
estimates around the critical point jsj]. In particular, for rigidity percolation on regular 
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lattices, the scenario of a second order phase transition in a different universality class than 
ordinary percolation seems favored by the numerics, although there has been some debate 
on the subject [10 j . 

From the analytical point of view however, progresses have been slow. A field theoretical 
attempt by Obukhov 11| predicts a first order phase transition in 2D, which, when com- 
pared with the numerics, seems to be a non generic feature. Some insight came from rigidity 
models exactly solved on trees and various types of random graphs, with locally tree-like 
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151 ] . However, in these types of models, the rigidity percolation tran- 



sition is also usually first order. This casts doubt on their usefulness to understand generic 
2D rigidity. There has been also earlier attempts to study rigidity percolation through 



Position Space Renormalization Group 
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171 ]; based on small renormalization cells, the 



associated predictions for the critical exponents are unprecise. To this date, we are not 
aware of any precise analytical prediction for the critical exponents of 2D generic rigidity 
percolation. 

A few years ago, Thorpe et al. have opened a new research direction by introducing 
the notion of network self-organization in rigidity percolation models . This is a natural 
idea in the context of network glasses modeling: constraints that are not satisfied create 
stress and bear an energy cost; the network should then self-organize, i.e. tend to modify 



its structure, in order to minimize this energy cost. In [18], numerical simulations allowing 
self organization point to the existence of a new phase in between the usual floppy and rigid 
ones. As a function of the mean connectivity, the phase diagram would then show two phase 
transitions instead of one. Around the same time, several ex per iments on network glasses 
drew a picture compati 



Die with this predicted phenomenology 



original simulations of 
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2l| . A problem of the 



181 ] was their strongly out of equilibrium character; some theoretical 
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23j. In |2J, [25|, a 



studies improved on this, and confirmed the "three phases" picture 
class of self-organizing rigidity percolation models is exactly solved, and shows three phases, 
separated by two true thermodynamic phase transitions. However, these exactly solvable 
models are based on random networks, which have a locally tree like topology, without small 
loops. This strongly influences the critical properties of the rigidity transition; in particular, 
as already mentioned above, rigidity percolation is continuous for a regular 2D triangular 



network, and it is first order in random networks. The simulations of refs. 
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271 ]. performed 



at zero temperature on a triangular lattice, seem to confirm the phenomenology suggested 
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by the locally tree like random networks studies, and show that the entire intermediate 



phase has a critical nature. Sartbaeva et al. in [28| discuss the influence of local structures 
on the existence and width of the intermediate phase in more realistic network glasses 
models. Despite these works, the question of whether self-organized rigidity percolation on 
two dimensional networks at finite temperature displays a true "intermediate phase" is still 
open. We conclude this paragraph by mentioning recent reviews on the subject (29), \^\ - 

The goal of this paper is twofold. First, we introduce a new class of solvable models 
of 2D rigidity, on hierarchical lattices. These lattices, in contrast with random graphs 
are not locally tree like, and possess many small loops. We may then expect for these 
lattices a phenomenology closer to generic 2D rigidity than random graphs. It is known 



that models defined on 
transformations 
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lierarchical lattices may be exactly solvable, using renormalization 



33 we will use the same ideas in the context of rigidity percolation. 



Second, in the spirit of 



24] for random networks, we add to the hierarchical networks 



the possibility of self-organization. This provides a new class of exactly solvable models, 
genuinely different from the random graphs of 24, |25|, where it is possible to investigate 
the existence of an "intermediate phase" . 

The paper is organized as follows. In section[TTl we summarize the combinatorial approach 
to 2D rigidity, which will be necessary in the following. After introducing, in section IHH 
the class of hierarchical lattices we will be interested in, we solve our model and study 
its critical properties in section IIVI We then introduce the possibility of self-organization 
for the hierarchical networks, and solve the associated model, in section IV] we discuss the 
implications on the existence of an "intermediate phase", between the usual floppy and 
stressed rigid ones. Section [VTl is devoted to discussion and perspectives. 



II. GENERIC 2D RIGIDITY 



An intuitive approach to rigidity percolation, apparently due to Maxwell, consists in 
counting degrees of freedom and constraints. In two dimensions, each site has two degrees 
of freedom; each link brings one constraint. In the end, a certain number of degrees of 
freedom are left. They correspond to unconstrained motions of the sites, and are usually 
called "floppy modes" . If one thinks of each link as a spring, with a given rest length, floppy 
modes are motions of the sites that do not change any spring length, and thus do not cost 
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any energy. Any network embedded in two dimensions necessarily has at least three floppy 
modes, corresponding to the two global translations and the global rotation in the plane. 
The network, or a subgraph of the network, is said to be rigid if it has no internal degree of 
freedom, besides these three. In a first approximation, a network with N sites and M links 
should have 

Nfiop = max(3, 2N - M) (1) 

degrees of freedom left, and be rigid if 2N — M < 3. This simple computation is obviously 
wrong in general, as shown on Fig. [TJ The reason is that when a new bond is added, it may 
be "redundant" : then the constraint associated with this new bond cannot be satisfied and 
its addition does not remove any degree of freedom. The exact formula giving the number 
of floppy modes of a network is 

N flop = 2N — M + N red , (2) 

where N re d is the number of redundant bonds. The problem is then reduced to the estimation 
of N ret i. The simplest form of Maxwell counting given by Eq. ([1]) assumes that new bonds 
added in the network first exhaust the degrees of freedom available before starting to become 
redundant. It is not correct, but it turns out that a generalization of this idea is actually 
exact 8|. It may be formulated as follows. Consider a graph, and add a new bond to it. It is 
redundant if and only if there exists a subgraph with N su b sites and M su b bonds, containing 
this new bond, such that 2N su b — M su b < 3. This combinatorial formulation is at the root of 
the "pebble game" , a powerful algorithm to analyze rigidity percolation on two dimensional 
networks 



III. HIERARCHICAL LATTICES 



We will consider the following type of hierarchical graphs. We start from two sites, 
connected by one bond. The graph is then constructed iteratively; at each step, all bonds 
are replaced by a given elementary cell. Four examples of elementary cells, corresponding to 
four examples of hierarchical graphs, are given on Fig. [21 In each case, it is easy to obtain a 
recursion relation for the number of sites N t and bonds B t of the graphs after t iterations. 
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FIG. 1: An example of a graph with N = 5 sites. Before the addition of the dashed bond, there 
are M = 6 bonds, none of which is redundant. The number of floppy modes is 4: two global 
translations, one global rotation, and the additional rotation of the one-coordinated site. The 
dashed bond is redundant: the subgraph isolated by the dashed curve has N BU f, = 4 sites and 
M su b = 6 > 2N su b — 3 bonds. The whole graph with the dashed bond also has 4 floppy modes, in 
agreement with Eq. (|2j). 



FIG. 2: Four examples of hierarchical graphs. On the left, the graphs at step t = 0, with two sites 
and one bond; on the right, the graphs at t = 1. We call the t = 1 graph the "elementary cell". 
1) is an elementary cell with 5 sites and 9 bonds. This is the hierarchical graph we will mostly 
consider in the following. 2), 3) and 4) are elementary cells with respectively 4 sites and 5 bonds, 
7 sites and 12 bonds, and 6 sites and 12 bonds. At step t, each bond of step t — 1 is replaced by a 
subgraph identical to the t = 1 graph. 

We obtain respectively, for the four graphs of Fig. [2] 
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Solving these recursion relations yields: 

3 v Q* 4- 1 3 

(7) 
(8) 
(9) 
(10) 

To model the rigidity percolation phenomenon, we now assume that each bond is effectively 
present in the graph with probability p; thus, the number of occupied bonds is, for large 
graphs, Mt — pB t . The Maxwell counting procedure yields an estimate of the threshold 
for rigidity percolation: the approximate critical probability puax corresponds to a ratio 
between bonds and sites M t /N t = 2. For graph on Fig. [212, we find p^ ax = 1; any p < 1 
should then render this graph floppy, and there is no true rigidity percolation. We thus 
dismiss this example in the following. For the graphs on Figs. [211, El3 and [2j4, we find 
respectively p$ ax = 3/4, p^ ax = 10/11 and p$ ax = 8/11. All three graphs should then 
present a rigidity transition at p < 1. We will mainly concentrate in the following on graph 
1) for the sake of simplicity. 



IV. SOLUTION OF THE HIERARCHICAL MODELS 

The goal is to compute a recursion relation relating N[ ed (p), the number of redundant 
bonds of a graph iterated t times, with a bond occupation probability p, and N™f(p'), 
the number of redundant bonds of a graph iterated t — 1 times, with a bond occupation 



probability p' . This will be possible because of the peculiar structure of the network [32|]. 

The procedure is somewhat reciprocal to the construction of the hierarchical network. 
Consider a graph iterated t times, with a bond probability p. Consider two sites A and 
B separated by one elementary cell, see Fig. [31 If they are rigidly connected through this 
cell, we replace the whole cell by a single bond. If they are not, we delete the whole cell. 
Repeating this procedure for all elementary cells, we construct a new graph, iterated t — 1 
times, with a bond occupation probability p' = tp(p). The function (p(p) is constructed by 
counting all configurations leading to a rigid connection between A and B. This process is 
detailed on Table HJ for the graph on Fig. [2jl. 




FIG. 3: Renormalization of an elementary cell, and counting of redundant bonds. Top left: sites 
A and B are rigidly connected, so that the elementary cell is renormalized into a bond; there is 
one redundant bond. Top right: sites A and B are rigidly connected, no redundant bond. Bottom 
left: sites A and B are not rigidly connected, no redundant bond. Bottom right: sites A and B 
are not rigidly connected, one redundant bond. 

Collecting the different contributions, we obtain 

tp{p) = p 9 + 9p 8 q + 30pV + 12pY + 3p 5 g 4 , (11) 

where q = 1 — p. We assume that the number of redundant bonds is an extensive quantity; 
neglecting subdominant contributions and probabilistic fluctuations, we then write 

Nr d (p) = N t r(p) , (12) 

and our goal is to compute r(p). Analyzing one renormalization step, we write 

Nr d ( P ) = N t ^(p)) + NUt,p) (13) 

where N' red {t,p) represents all the redundant bonds that where suppressed in the renor- 
malization process. This formula needs some justification. Consider the empty (without 
bonds) graph at step t, and add the bonds one by one, checking each time if the newly 
added bond is redundant or not; this is actually the numerical strategy implemented in the 
"pebble game" algorithm. If the added bond is redundant, it is discarded. Thus, during 
the process, the constructed graph has at most one redundant bond. If an added bond is 
redundant, it is possible to find a subgraph S, containing this new bond, with N su b sites and 
M su b = 2N su b — 2 bonds. We choose S as small as possible. There are now two possibilities. 
Either S is included in one elementary cell; in this case the redundant bond is included in 
the N' red {t,p) contribution. Or S is not included in an elementary cell. It is easy to see that 
no floppy elementary cell can be included in S, so that the newly added bond necessarily 
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TABLE I: Enumeration of all types of possible configurations for an elementary cell of the graph 
on Fig. [2jl. There are for instance 12 configurations with 6 bonds, providing a rigid connection 
across the cell, and without redundant bond. Some examples are given on Fig. 

"rigidities" the elementary cell it belongs to. In addition, if S contains a bond included in 
one elementary cells, S contains both sites which are end points of the cell, and the cell is 
rigid. S, which is included in the graph at step t will then be renormalized in a natural 
way in a subgraph of the graph at step t — 1. The redundant constraint is thus included in 
the N™f(ip(p)) contribution to Eq. (jl~3]) . It is worth mentioning that both the hierarchical 
nature of the graph and the convenient combinatorial description of 2D rigidity outlined in 
Section [TT] are necessary for this step. In particular, it is not clear how to construct exactly 
solvable hierarchical models of 3D rigidity. 

For a given bond probability p, we want to compute r (p), the mean number of redundant 
bonds inside one elementary cell. We use Table [T]to enumerate all possibilities, and we obtain 
the following expression for ro(p): 

r (p) = 2p 9 + 9p s q + 6p 7 q 2 + 2p 6 q 3 (14) 
9 




FIG. 4: (Color online) The red curve y = <p(p) (thick line), together with the blue y = p straight 
line (thin line). The critical point defined by <p(p*) = p* is indicated by an arrow. 

In one renormalization step, a large number of elementary cells B t -\ is renormalized. 
Thus N' red {t,p) ~ B t -iro(p). Finally, the recursion relation reads 



N t r(p) = N t -ir(<p(p)) + B t _ir (p) (15) 

Using the relations B t = B l and B t /N t = x , and simplifying by N t , this yields 

r(p) = — r (ip(p)) + f (p) , (16) 

with f (p) = xqTq(p)/ B. For our main example, represented on Fig. [3l B = 9 and Xq = 8/3. 

From Eq. (121), we see that for given numbers of sites and bonds, the number of flo ppy 
modes is linearly related to the number of redundant bonds. Since it has been noticed 10] 
that the number of floppy modes plays the role of a "free energy" in rigidity percolation 
models, we remark that the computation of r(p) is analogous to the free energy computation 



of ref. [33j, for spin models on hierarchical lattices. 

The function if is plotted on Fig. |H It has three fixed points; p = and p — 1 are stable 
and correspond to the floppy and rigid phases respectively, p = p* ~ 0.867 is unstable 
and is the critical point. We note ip^ the fcth iterate of <p. Then for all p < p* (resp. 
p > p*), Lp( n \p) tends to (resp. 1). The quantity a = dip/dp(p = p*) ~ 1.97 will play an 
important role. We note that if we had used the hierarchical network on Fig. [2j2, we would 
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have obtained a function tp with only two fixed points: p = 0, stable, and p = 1, unstable. 
Thus, in this case, any p < 1 yields a floppy network at large scales. A similar analysis for 
the graphs on Figs. [213 and [214 yields respectively p* = 0.975, p* = 0.853 for the critical 
densities, and a = 1.86, a = 2.07 for the tp derivative at criticality. We note that in all three 
cases, we find p* > puax- 

Iterating n times Eq. (|T6|) . one obtains 



rip) = ^r^\p)) + ^- T r ^- 1 \p)) + ...+ 

+ -g f ((p(p)) + f (p) (17) 

As B > 1 and f is bounded, r may then be expressed as an infinite series: 

oo 1 
fc=0 

As the k times iterated function (p^ k '(p) tends pointwise to when p < p*, and to 1 when 
p > p*, the above series is infinitely differentiable, except when p = p*, the unstable fixed 
point of tp. Looking for a power law behavior of r close to p* 

r(p) = Cte\p — p*\ T + regular part , (19) 

one finds r = In Bj In a. For our main example, we find r = 3.24; this corresponds to 
a critical exponent a = —1.24, where we have used the standard notation for the critical 
exponent associated with the free energy singularity. This means that the derivatives of r(p) 
are all continuous up to third order, while the fourth derivative is discontinuous. This does 
not compare well to the numerical results on triangular networks 9J, which find a ~ —0.48. 
Networks on Fig.[2l3 and Fig.[2l4 yield respectively r = 4.01 and r = 3.41, which corresponds 
to a = —2.01 and a = —1.41. Thus, there are important fluctuations in the value of the 
critical exponent between two different hierarchical networks. We note however that the 
tendency to a very weak rigidity percolation transition on these hierarchical networks seems 
to be general. A natural idea to improve the agreement between these analytical results on 
hierarchical networks and the numerics on triangular lattices would be to consider larger 
renormalization cells. The previous remarks cast doubt on this strategy. Another interesting 
extension would be to study other critical properties, like the size of the largest rigid cluster, 
or the largest stressed cluster. 
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V. ADAPTIVE NETWORKS 



A. Introduction of an energy 

We consider now the possibility for the network to modify its structure, in order to avoid 
stress. This step was first taken in [l8], in an attempt to better model network glasses; in 
this context, it is natural to assume that the connectivity structure of the network may adapt 
itself to a certain extent to avoid stressed bonds. Following , we introduce an energy cost 
for each constraint which is not satisfied; at any finite temperature, this cost competes with 
the entropic cost of reorganizing the network. In this setting, standard rigidity percolation, 
as studied up to now, corresponds to the infinite temperature limit, where the energy cost 
is negligible. 

We have to define an energy for any configuration of the network. The simplest choice 
is to count the number of redundant bonds, and to attribute the same amount of energy to 



each of them 



24 



25J. 



We consider a network where a proportion x of the bonds is present. The partition 
function for our model, on a hierarchical network iterated t times, reads: 

Z t (x,P)= e-^W'-H (20) 

{Ji},E*i=*-Bt 

where we note lj = 1 if bond j is present, lj = otherwise; Yl{i-y * s the sum over all possible 
configurations of the network. Here, it is restricted to the configurations with the prescribed 
fraction of occupied bonds, x. 

To proceed, it is useful to relax the constraint on the number of bonds. We introduce a 
"grand canonical" statistical ensemble, whose partition sum reads: 

where p is the a priori probability that a bond is present, and q = 1 — p. ^2 lj is the total 
number of bonds in the network in configuration {lj}. We recall that B t is the maximum 
number of bonds. Because of the energetic term which biases the probability distribution 
towards non stressed configurations, the actual number of bonds present in the network is 
different from pB t , except in the infinite temperature limit. We may rewrite the partition 
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sum as 



n t (p >/3 ) = (j^j '^e-^e-^WV, (22) 



where 



A = In ( ^—^ ) (23) 



P 

We define ou(p, P) = hiQ t {p, f3)/N t . Generalizing the renormalization transformation of the 
previous section, it is possible to compute exactly u(p,(3). We consider one elementary cell, 
and want to compute <p{p,j3), the probability that this elementary cell is renormalized as 
a bond. The calculation is the same as in section IIVI except that we have to attribute a 
weight e~ n P instead of 1 to configurations with n redundant bonds. For the elementary cell 
of our chief example, Fig. [2jl, we use once again Table [U to obtain 

<pi(p, P) + <Po[p, P) 

where tpi collects the contributions of the configurations where the end sites of the cell are 
rigidly connected, and (fo collects the contributions of the configurations where the end sites 
of the cell are not rigidly connected. Division by <pi + tfo is for normalization. We have 

^(p, p) = p 9 e ~ W + 9p 8 qe~ 13 + 30p 7 q 2 + 12p 6 g 3 + 3p 5 q 4 
^ ( p? p) = Qp 7 q 2 e ~ f3 + 2p 6 q 3 e-P + 70p 6 q 3 + 123pV 

+126pV + 84p 3 g 6 + 36p 2 q 7 + 9pq 8 + q 9 (25) 

Some curves <p(p, P) are drawn for various values of P on Fig. We also define 

c(p,P)=tpi(p t P)+tp (p,P) (26) 

Then, following the same reasoning as above, we derive the recursion relation 

Qt(p,P) = c(p,/3) B *- 1 ^_ 1 (^(p,/3),/3) (27) 

Taking the logarithm: 

N t u(p,P) = B t . 1 hic(p,P) + N t -. 1 u;(<P(p,P),P) ( 2 §) 
Recalling that N t — 3 X 9*/8 and B t = 9', we have 

u(p,P) = ^kic(p,P) + ^u(ip(p,P),P) (29) 
13 
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FIG. 5: (Color online) Examples of functions tp(p,/3), as a function of p; from left to right, (3 = 
0, 1, 5 (thick lines, respectively red, black and light blue curves). The intersection with the y = p 
straight line (blue thin line) is the critical point. A low temperature shifts the critical point towards 
high p values. 

This yields the following expression for uj 



As above, ip( k > is the kth. iterate of the function </?. We note that the recursion relation 
for p depends parametrically on j3, and that the temperature is not renormalized in the 
process. This implies that there exists a line of critical points, and that the critical exponents 
continuously depend on (3 along this line. From Eqs. fl22l) and fl23|) . we obtain an expression 
for x(p, (3) 



We would like to study our model at fixed x; the previous formula tells us how to choose 
p to do so, providing the bridge between canonical and grand canonical ensembles. Fig. [6] 
shows various x(p) curves, for different values of /3. From Fig. [5], we can see that p*(/3) tends 
to 1 for large j3, as we could have anticipated. However, the critical connectivity x*(/3) does 
not tend to 8/3, see Fig. [61 A crucial observable is the energy per site, that is the number 




(30) 



duo 
dp 



(31) 
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FIG. 6: (Color online) The mean connectivity of the network x as a function of p. From left to 
right, (3 = 0, 1, 5 (respectively red, black and light blue curves). The critical points are also 
indicated. 

of redundant bonds per site n re d] it is computed using the formula: 

dui . . 

n r ed = (32) 

n red(p, P) curves are shown on Fig. [TJ together with results of Monte Carlo simulations. 



B. An intermediate phase? 

From the exact solution of the model provided in the previous paragraph, it is clear 
that the introduction of an energy and a finite temperature does not induce any qualitative 
change in the phase diagram: as a function of the connectivity, there are only two phases, 
floppy and rigid, separated by a phase transition at a critical connectivity x*(/3). Thus, 
our first conclusion is that in this model, the "intermediate phase" does not have a real 
thermodynamical meaning. It is an important difference with the case of the random net- 



works 



24 



25j . This result seems related to the strength of the percolation transition, which 
is first order in random networks, and very weak on the hierarchical lattices introduced in 
this paper. 

However, looking at the curves on Fig. 0, it is possible at low temperatures to distinguish 

15 
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FIG. 7: (Color online) The number of redundant per site is plotted as a function of the mean 
connectivity x, for (3 = 0, 1, 5 (from left to right, respectively red, black and light blue curves). 
The crosses indicate the critical points; the filled circles are results of Metropolis Monte Carlo 
simulations, for (3 = 1, using the pebble game algorithm, on hierarchical graphs with t = 5 
(22145 sites). They are in perfect agreement with the analytical results. 

three regions for low enough temperatures. At low connectivity x < 2, there is essentially no 
redundant bond, and thus no stress in the network; for 2 < x < x*((3), there is some stress 
in the network, but it does not percolate; finally, for x > x*(/3), the stress percolates in the 
network. The boundary between the two first regions is not a ther mo dynamical transition; 
rather, it is a cross-over which becomes sharper as the temperature decreases. In this cross- 
over region, isostatic local structures are favored, and the effect of network self-organization 
is the strongest. To illustrate these points, we have plotted on Fig. [8] the fraction of isostatic 
elementary cells (with exactly 7 bonds). To define an indicator for self-organization, we 
note that without network adaptivity (or equivalently at infinite temperature), the mean 
connectivity x is directly proportional to the parameter p, and given by x = 8p/3; we may 
then define a self-organization indicator as s = (p—3x/8)/p. s is plotted on Fig. [8], and peaks 
in the cross over region, just as the fraction of elementary cells. The lower the temperature, 
the sharper the peaks get. 

Finally, we stress that we find a percolation transition at a critical density x c ~ 2.3, 
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FIG. 8: (Color online) Top panel: the self organization indicator s = (p — 3x/8)/p, as a function of 
the mean connectivity x. Bottom panel: fraction of isostatic elementary cells as a function of the 
mean connectivity. In both panels, (3 = 0, 1, 5 from bottom to top (respectively red, black and 
light blue curves), and the crosses indicate the critical points. Note that s is identically when 
(3 = 0. 

larger than the Maxwell threshold XMax = 2. This contrasts with rigidity models on random 
graphs, for which the critical density is usually lower than xuax- 
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VI. CONCLUSION 



Our first result is the construction of solvable models for rigidity percolation, beyond 
trees and random networks. These models display a critical point at the rigidity percolation 
threshold, and the critical exponent for the number of redundant bonds (which corresponds 
to the free energy) can be computed exactly. It would be interesting to compute also the 
exponent associated with the biggest rigid cluster in the network. In any case, it is known 
that these exponents, computed on such hierarchical networks, cannot be taken as reliable 
estimates for the exponents associated to regular lattices. However, these exponents may 
converge to the correct ones when one considers larger and larger cells for the hierarchical 



networks 



34|. One may then think of a large cell Monte Carlo renormalization procedure 



to obtain an estimation of the critical exponents for 2D generic rigidity percolation. Some 
words of caution are in order: it is not clear how to extend to large sizes the hierarchical 
networks considered in this paper, nor if they may be a good approximation of a regular 
lattice, even for large cell sizes. It may be necessary to introduce other types of hierarchical 
lattices. 

Our second result concerns the effect of self organization on rigidity percolation, for these 
hierarchical models. We define and solve exactly a model where the network may adapt its 
structure to avoid stress, at an entropic cost. At variance with what happens for random 
networks, the possibility of self-organization does not introduce a qualitative change in the 
phase diagram. In particular, there is no true thermodynamic "intermediate phase" . It may 
still be possible to define an "intermediate region" , where the network is close to be isostatic, 
and the stress does not percolate, although it is present locally. It would be interesting to 

re picture of the intermediate 



22j. From the modeling point 



compare the properties of this "intermediate region" with t 
phase drawn from a Size Increasing Cluster Approximation 
of view, we note that a scenario where sharp phase transitions are replaced by smoother 
cross-overs may still be compatible with the experiments on network glasses. 

All our analysis relies on the combinatorial description of 2D generic rigidity. It should be 



possible to study similar hierarchical models o: 
using the associated combinatorial description 



3D rigidity with bond bending constraints, 



35[ . It is not clear however how to generalize 



the technique to 3D networks with central forces only, as there is no simple combinatorial 
description in that case (see 36( for a recent discussion). 
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